A pilot investigation of differential hydroxymethylation levels in patient-derived neural stem cells implicates altered cortical development in bipolar disorder

IntroductionBipolar disorder (BD) is a chronic mental illness characterized by recurrent episodes of mania and depression and associated with social and cognitive disturbances. Environmental factors, such as maternal smoking and childhood trauma, are believed to modulate risk genotypes and contribute to the pathogenesis of BD, suggesting a key role in epigenetic regulation during neurodevelopment. 5-hydroxymethylcytosine (5hmC) is an epigenetic variant of particular interest, as it is highly expressed in the brain and is implicated in neurodevelopment, and psychiatric and neurological disorders.MethodsInduced pluripotent stem cells (iPSCs) were generated from the white blood cells of two adolescent patients with bipolar disorder and their same-sex age-matched unaffected siblings (n = 4). Further, iPSCs were differentiated into neuronal stem cells (NSCs) and characterized for purity using immuno-fluorescence. We used reduced representation hydroxymethylation profiling (RRHP) to perform genome-wide 5hmC profiling of iPSCs and NSCs, to model 5hmC changes during neuronal differentiation and assess their impact on BD risk. Functional annotation and enrichment testing of genes harboring differentiated 5hmC loci were performed with the online tool DAVID.ResultsApproximately 2 million sites were mapped and quantified, with the majority (68.8%) located in genic regions, with elevated 5hmC levels per site observed for 3’ UTRs, exons, and 2-kb shorelines of CpG islands. Paired t-tests of normalized 5hmC counts between iPSC and NSC cell lines revealed global hypo-hydroxymethylation in NSCs and enrichment of differentially hydroxymethylated sites within genes associated with plasma membrane (FDR = 9.1 × 10−12) and axon guidance (FDR = 2.1 × 10−6), among other neuronal processes. The most significant difference was observed for a transcription factor binding site for the KCNK9 gene (p = 8.8 × 10−6), encoding a potassium channel protein involved in neuronal activity and migration. Protein–protein-interaction (PPI) networking showed significant connectivity (p = 3.2 × 10−10) between proteins encoded by genes harboring highly differentiated 5hmC sites, with genes involved in axon guidance and ion transmembrane transport forming distinct sub-clusters. Comparison of NSCs of BD cases and unaffected siblings revealed additional patterns of differentiation in hydroxymethylation levels, including sites in genes with functions related to synapse formation and regulation, such as CUX2 (p = 2.4 × 10−5) and DOK-7 (p = 3.6 × 10−3), as well as an enrichment of genes involved in the extracellular matrix (FDR = 1.0 × 10−8).DiscussionTogether, these preliminary results lend evidence toward a potential role for 5hmC in both early neuronal differentiation and BD risk, with validation and more comprehensive characterization to be achieved through follow-up study

A pilot investigation of differential hydroxymethylation levels in patient-derived neural stem cells implicates altered cortical development in bipolar disorder

Introduction: Bipolar disorder (BD) is a chronic mental illness characterized by recurrent episodes of mania and depression and associated with social and cognitive disturbances. Environmental factors, such as maternal smoking and childhood trauma, are believed to modulate risk genotypes and contribute to the pathogenesis of BD, suggesting a key role in epigenetic regulation during neurodevelopment. 5-hydroxymethylcytosine (5hmC) is an epigenetic variant of particular interest, as it is highly expressed in the brain and is implicated in neurodevelopment, and psychiatric and neurological disorders. Methods: Induced pluripotent stem cells (iPSCs) were generated from the white blood cells of two adolescent patients with bipolar disorder and their same-sex age-matched unaffected siblings (n = 4). Further, iPSCs were differentiated into neuronal stem cells (NSCs) and characterized for purity using immuno-fluorescence. We used reduced representation hydroxymethylation profiling (RRHP) to perform genome-wide 5hmC profiling of iPSCs and NSCs, to model 5hmC changes during neuronal differentiation and assess their impact on BD risk. Functional annotation and enrichment testing of genes harboring differentiated 5hmC loci were performed with the online tool DAVID. Results: Approximately 2 million sites were mapped and quantified, with the majority (68.8%) located in genic regions, with elevated 5hmC levels per site observed for 3\u27 UTRs, exons, and 2-kb shorelines of CpG islands. Paired t-tests of normalized 5hmC counts between iPSC and NSC cell lines revealed global hypo-hydroxymethylation in NSCs and enrichment of differentially hydroxymethylated sites within genes associated with plasma membrane (FDR = 9.1 × 10-12) and axon guidance (FDR = 2.1 × 10-6), among other neuronal processes. The most significant difference was observed for a transcription factor binding site for the KCNK9 gene (p = 8.8 × 10-6), encoding a potassium channel protein involved in neuronal activity and migration. Protein-protein-interaction (PPI) networking showed significant connectivity (p = 3.2 × 10-10) between proteins encoded by genes harboring highly differentiated 5hmC sites, with genes involved in axon guidance and ion transmembrane transport forming distinct sub-clusters. Comparison of NSCs of BD cases and unaffected siblings revealed additional patterns of differentiation in hydroxymethylation levels, including sites in genes with functions related to synapse formation and regulation, such as CUX2 (p = 2.4 × 10-5) and DOK-7 (p = 3.6 × 10-3), as well as an enrichment of genes involved in the extracellular matrix (FDR = 1.0 × 10-8). Discussion: Together, these preliminary results lend evidence toward a potential role for 5hmC in both early neuronal differentiation and BD risk, with validation and more comprehensive characterization to be achieved through follow-up study

Whole Genome Sequence Data From Captive Baboons Implicate RBFOX1 in Epileptic Seizure Risk

In this study, we investigate the genetic determinants that underlie epilepsy in a captive baboon pedigree and evaluate the potential suitability of this non-human primate model for understanding the genetic etiology of human epilepsy. Archived whole-genome sequence data were analyzed using both a candidate gene approach that targeted variants in baboon homologs of 19 genes (n = 20,881 SNPs) previously implicated in genetic generalized epilepsy (GGE) and a more agnostic approach that examined protein-altering mutations genome-wide as assessed by snpEff (n = 36,169). Measured genotype association tests for baboon cases of epileptic seizure were performed using SOLAR, as well as gene set enrichment analyses (GSEA) and protein–protein interaction (PPI) network construction of top association hits genome-wide (p \u3c 0.01; n = 441 genes). The maximum likelihood estimate of heritability for epileptic seizure in the pedigreed baboon sample is 0.76 (SE = 0.77; p = 0.07). Among candidate genes for GGE, a significant association was detected for an intronic SNP in RBFOX1 (p = 5.92 × 10–6; adjusted p = 0.016). For protein-altering variants, no genome-wide significant results were observed for epilepsy status. However, GSEA revealed significant positive enrichment for genes involved in the extracellular matrix structure (ECM; FDR = 0.0072) and collagen formation (FDR = 0.017), which was reflected in a major PPI network cluster. This preliminary study highlights the potential role of RBFOX1 in the epileptic baboon, a protein involved in transcriptomic regulation of multiple epilepsy candidate genes in humans and itself previously implicated in human epilepsy, both focal and generalized. Moreover, protein-damaging variants from across the genome exhibit a pattern of association that links collagen-containing ECM to epilepsy risk. These findings suggest a shared genetic etiology between baboon and human forms of GGE and lay the foundation for follow-up research

Dopamine perturbation of gene co-expression networks reveals differential response in schizophrenia for translational machinery.

The dopaminergic hypothesis of schizophrenia (SZ) postulates that positive symptoms of SZ, in particular psychosis, are due to disturbed neurotransmission via the dopamine (DA) receptor D2 (DRD2). However, DA is a reactive molecule that yields various oxidative species, and thus has important non-receptor-mediated effects, with empirical evidence of cellular toxicity and neurodegeneration. Here we examine non-receptor-mediated effects of DA on gene co-expression networks and its potential role in SZ pathology. Transcriptomic profiles were measured by RNA-seq in B-cell transformed lymphoblastoid cell lines from 514 SZ cases and 690 controls, both before and after exposure to DA ex vivo (100 μM). Gene co-expression modules were identified using Weighted Gene Co-expression Network Analysis for both baseline and DA-stimulated conditions, with each module characterized for biological function and tested for association with SZ status and SNPs from a genome-wide panel. We identified seven co-expression modules under baseline, of which six were preserved in DA-stimulated data. One module shows significantly increased association with SZ after DA perturbation (baseline: P = 0.023; DA-stimulated: P = 7.8 × 10-5; ΔAIC = -10.5) and is highly enriched for genes related to ribosomal proteins and translation (FDR = 4 × 10-141), mitochondrial oxidative phosphorylation, and neurodegeneration. SNP association testing revealed tentative QTLs underlying module co-expression, notably at FASTKD2 (top P = 2.8 × 10-6), a gene involved in mitochondrial translation. These results substantiate the role of translational machinery in SZ pathogenesis, providing insights into a possible dopaminergic mechanism disrupting mitochondrial function, and demonstrates the utility of disease-relevant functional perturbation in the study of complex genetic etiologies

Variants in the human potassium channel gene (KCNN3) are associated with migraine in a high risk genetic isolate

The calcium-activated potassium ion channel gene (KCNN3) is located in the vicinity of the familial hemiplegic migraine type 2 locus on chromosome 1q21.3. This gene is expressed in the central nervous system and plays a role in neural excitability. Previous association studies have provided some, although not conclusive, evidence for involvement of this gene in migraine susceptibility. To elucidate KCNN3 involvement in migraine, we performed gene-wide SNP genotyping in a high-risk genetic isolate from Norfolk Island, a population descended from a small number of eighteenth century Isle of Man ‘Bounty Mutineer’ and Tahitian founders. Phenotype information was available for 377 individuals who are related through the single, well-defined Norfolk pedigree (96 were affected: 64 MA, 32 MO). A total of 85 SNPs spanning the KCNN3 gene were genotyped in a sub-sample of 285 related individuals (76 affected), all core members of the extensive Norfolk Island ‘Bounty Mutineer’ genealogy. All genotyping was performed using the Illumina BeadArray platform. The analysis was performed using the statistical program SOLAR v4.0.6 assuming an additive model of allelic effect adjusted for the effects of age and sex. Haplotype analysis was undertaken using the program HAPLOVIEW v4.0. A total of four intronic SNPs in the KCNN3 gene displayed significant association (P < 0.05) with migraine. Two SNPs, rs73532286 and rs6426929, separated by approximately 0.1 kb, displayed complete LD (r2 = 1.00, D′ = 1.00, D′ 95% CI = 0.96–1.00). In all cases, the minor allele led to a decrease in migraine risk (beta coefficient = 0.286–0.315), suggesting that common gene variants confer an increased risk of migraine in the Norfolk pedigree. This effect may be explained by founder effect in this genetic isolate. This study provides evidence for association of variants in the KCNN3 ion channel gene with migraine susceptibility in the Norfolk genetic isolate with the rarer allelic variants conferring a possible protective role. This the first comprehensive analysis of this potential candidate gene in migraine and also the first study that has utilised the unique Norfolk Island large pedigree isolate to implicate a specific migraine gene. Studies of additional variants in KCNN3 in the Norfolk pedigree are now required (e.g. polyglutamine variants) and further analyses in other population data sets are required to clarify the association of the KCNN3 gene and migraine risk in the general outbred population

Methylation of SOCS3 is inversely associated with metabolic syndrome in an epigenome-wide association study of obesity

Epigenetic mechanisms, including DNA methylation, mediate the interaction between gene and environment and may play an important role in the obesity epidemic. We assessed the relationship between DNA methylation and obesity in peripheral blood mononuclear cells (PBMCs) at 485,000 CpG sites across the genome in family members (8-90 y of age) using a discovery cohort (192 individuals) and a validation cohort (1,052 individuals) of Northern European ancestry. After Bonferroni-correction (Pα=0.05 = 1.31 × 10-7) for genome-wide significance, we identified 3 loci, cg18181703 (SOCS3), cg04502490 (ZNF771), and cg02988947 (LIMD2), where methylation status was associated with body mass index percentile (BMI%), a clinical index for obesity in children, adolescents, and adults. These sites were also associated with multiple metabolic syndrome (MetS) traits, including central obesity, fat depots, insulin responsiveness, and plasma lipids. The SOCS3 methylation locus was also associated with the clinical definition of MetS. In the validation cohort, SOCS3 methylation status was found to be inversely associated with BMI% (P = 1.75 × 10-6), waist to height ratio (P = 4.18 × 10-7), triglycerides (P = 4.01 × 10-4), and MetS (P = 4.01 × 10-7), and positively correlated with HDL-c (P = 4.57 × 10-8). Functional analysis in a sub cohort (333 individuals) demonstrated SOCS3 methylation and gene expression in PBMCs were inversely correlated (P = 2.93 × 10-4) and expression of SOCS3 was positively correlated with status of MetS (P = 0.012). We conclude that epigenetic modulation of SOCS3, a gene involved in leptin and insulin signaling, may play an important role in obesity and MetS

Whole genome sequence data implicate RBFOX1 in epilepsy risk in baboons

Background: Baboons exhibit a genetic generalized epilepsy (GGE) that resembles juvenile myoclonic epilepsy and may represent a suitable genetic model for human epilepsy. The genetic underpinnings of epilepsy were investigated in a baboon colony at the Southwest National Primate Research Center (San Antonio, TX) through the analysis of whole-genome sequence (WGS) data. Methods: Baboon WGS data were obtained for 38 cases and 19 healthy controls from the NCBI Sequence Read Archive and, after standard QC filtering, two subsets of variants were examined: (1) 20,881 SNPs from baboon homologs of 19 candidate GGE genes; and (2) 36,169 protein-altering SNPs. Association tests were conducted in SOLAR, and gene set enrichment analyses (GSEA) and protein-protein interaction (PPI) network construction were performed on genome-wide significant association results (Pn= 441 genes). Results: Heritability for epileptic seizure in the pedigreed baboon sample was estimated at 0.76 (SE=0.77; P=0.07). A significant association was detected for an intronic SNP in RBFOX1 (P=5.92 × 10-6; adjusted P=0.016). For protein-altering variants, GSEA revealed significant positive enrichment for genes involved in the extracellular matrix structure (ECM; FDR=0.0072) and collagen formation (FDR=0.017). Conclusions: SNP association results implicate RBFOX1 in baboon epilepsy, a gene that plays a key role in neuronal excitation and transcriptomic regulation, and has been previously linked to human epilepsy, both focal and generalized. Moreover, protein-damaging variants from across the baboon genome exhibit a wider pattern of association that links collagen-containing ECM to epilepsy risk. These findings suggest a shared genetic etiology between baboon and human forms of GGE

Association of differential gene expression with imatinib mesylate and omacetaxine mepesuccinate toxicity in lymphoblastoid cell lines

BackgroundImatinib mesylate is currently the drug of choice to treat chronic myeloid leukemia. However, patient resistance and cytotoxicity make secondary lines of treatment, such as omacetaxine mepesuccinate, a necessity. Given that drug cytotoxicity represents a major problem during treatment, it is essential to understand the biological pathways affected to better predict poor drug response and prioritize a treatment regime.MethodsWe conducted cell viability and gene expression assays to determine heritability and gene expression changes associated with imatinib and omacetaxine treatment of 55 non-cancerous lymphoblastoid cell lines, derived from 17 pedigrees. In total, 48,803 transcripts derived from Illumina Human WG-6 BeadChips were analyzed for each sample using SOLAR, whilst correcting for kinship structure.ResultsCytotoxicity within cell lines was highly heritable following imatinib treatment (h2&thinsp;=&thinsp;0.60-0.73), but not omacetaxine treatment. Cell lines treated with an IC20 dose of imatinib or omacetaxine showed differential gene expression for 956 (1.96%) and 3,892 transcripts (7.97%), respectively; 395 of these (0.8%) were significantly influenced by both imatinib and omacetaxine treatment. k-means clustering and DAVID functional annotation showed expression changes in genes related to kinase binding and vacuole-related functions following imatinib treatment, whilst expression changes in genes related to cell division and apoptosis were evident following treatment with omacetaxine. The enrichment scores for these ontologies were very high (mostly &gt;10).ConclusionsInduction of gene expression changes related to different pathways following imatinib and omacetaxine treatment suggests that the cytotoxicity of such drugs may be differentially tolerated by individuals based on their genetic background.<br /

Blood-Based miRNA Biomarkers as Correlates of Brain-Based miRNA Expression

The use of easily accessible peripheral samples, such as blood or saliva, to investigate neurological and neuropsychiatric disorders is well-established in genetic and epigenetic research, but the pathological implications of such biomarkers are not easily discerned. To better understand the relationship between peripheral blood- and brain-based epigenetic activity, we conducted a pilot study on captive baboons (Papio hamadryas) to investigate correlations between miRNA expression in peripheral blood mononuclear cells (PBMCs) and 14 different cortical and subcortical brain regions, represented by two study groups comprised of 4 and 6 animals. Using next-generation sequencing, we identified 362 miRNAs expressed at ≥ 10 read counts in 80% or more of the brain samples analyzed. Nominally significant pairwise correlations (one-sided P \u3c 0.05) between peripheral blood and mean brain expression levels of individual miRNAs were observed for 39 and 44 miRNAs in each group. When miRNA expression levels were averaged for tissue type across animals within the groups, Spearman\u27s rank correlations between PBMCs and the brain regions are all highly significant (r s = 0.47-0.57; P \u3c 2.2 × 10-16), although pairwise correlations among the brain regions are markedly stronger (r s = 0.86-0.99). Principal component analysis revealed differentiation in miRNA expression between peripheral blood and the brain regions for the first component (accounting for ∼75% of variance). Linear mixed effects modeling attributed most of the variance in expression to differences between miRNAs (\u3e70%), with non-significant 7.5% and 13.1% assigned to differences between blood and brain-based samples in the two study groups. Hierarchical UPGMA clustering revealed a major co-expression branch in both study groups, comprised of miRNAs globally upregulated in blood relative to the brain samples, exhibiting an enrichment of miRNAs expressed in immune cells (CD14+, CD15+, CD19+, CD3+, and CD56 + leukocytes) among the top blood-brain correlates, with the gene MYC, encoding a master transcription factor that regulates angiogenesis and neural stem cell activation, representing the most prevalent miRNA target. Although some differentiation was observed between tissue types, these preliminary findings reveal wider correlated patterns between blood- and brain-expressed miRNAs, suggesting the potential utility of blood-based miRNA profiling for investigating by proxy certain miRNA activity in the brain, with implications for neuroinflammatory and c-Myc-mediated processes

Dopamine perturbation of gene co-expression networks reveals differential response in schizophrenia for translational machinery

The dopaminergic hypothesis of schizophrenia (SZ) postulates that positive symptoms of SZ, in particular psychosis, are due to disturbed neurotransmission via the dopamine (DA) receptor D2 (DRD2). However, DA is a reactive molecule that yields various oxidative species, and thus has important non-receptor-mediated effects, with empirical evidence of cellular toxicity and neurodegeneration. Here we examine non-receptor-mediated effects of DA on gene co-expression networks and its potential role in SZ pathology. Transcriptomic profiles were measured by RNA-seq in B-cell transformed lymphoblastoid cell lines from 514 SZ cases and 690 controls, both before and after exposure to DA ex vivo (100 μM). Gene co-expression modules were identified using Weighted Gene Co-expression Network Analysis for both baseline and DA-stimulated conditions, with each module characterized for biological function and tested for association with SZ status and SNPs from a genome-wide panel. We identified seven co-expression modules under baseline, of which six were preserved in DA-stimulated data. One module shows significantly increased association with SZ after DA perturbation (baseline: P = 0.023; DA-stimulated: P = 7.8 × 10-5; ΔAIC = −10.5) and is highly enriched for genes related to ribosomal proteins and translation (FDR = 4 × 10−141), mitochondrial oxidative phosphorylation, and neurodegeneration. SNP association testing revealed tentative QTLs underlying module co-expression, notably at FASTKD2 (top P = 2.8 × 10−6), a gene involved in mitochondrial translation. These results substantiate the role of translational machinery in SZ pathogenesis, providing insights into a possible dopaminergic mechanism disrupting mitochondrial function, and demonstrates the utility of disease-relevant functional perturbation in the study of complex genetic etiologies